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We explain in this paper how a meaningful irrelevant perturbation theory around the 
infra-red (strong coupling) fixed point can be carried out for integrable quantum impurity 
problems. This is illustrated in details for the spin 1/2 Kondo model, where our approach 
gives rise to the complete low temperature expansion of the resistivity, beyond the well 
known T 2 Fermi liquid behaviour. We also consider the edge states tunneling problem, and 
demonstrate by Keldysh techniques that the DC current satisfies an exact duality between 
the UV and IR regimes. This corresponds physically to a duality between the tunneling 
of Laughlin quasi particles and electrons, and, more formally, to the existence of an exact 
instantons expansion. The duality is deeply connected with integrability, and could not 
have been expected a priori. 
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1. Introduction 

Duality arguments have been commonly used in quantum impurity problems for many 
years. An archetypal situation is provided by the model of a particle moving in a periodic 
potential and subject to quantum dissipation 0]. This problem is represented in one 
dimension @ by the following (boundary sine-Gordon or BSG) hamiltonian 

-i rO poo 

H = - dx dy[(d x $) 2 + (n) 2 ]+2\cosy/2^j$(0). (1.1) 

^ J — oo J — oo 

where the field at the origin represents the particle coordinate, and the bulk free boson 
the bath degrees of freedom. At small A the particle diffuses freely (UV fixed point) , 
while at large A, it is localized (IR fixed point) in a minimum of the potential, given 



by $(0) = \l g n i n an integer. Near the UV fixed point, physical properties can be 
expanded in powers of A, and are expressed in terms of Coulomb gas integrals whose 
charges ±g correspond to the two possible exponentials in the cosine term (the dimension 



of the operator cos y/2ng&(0) being A = g). It is also possible to study the vicinity of 
this IR fixed point in an \ expansion by considering the instantons and anti-instantons 
that take the particle from one minima to a neighbouring one. Using the leading order 
action of these instantons, one obtains again a Coulomb gas, but this time with charges 
±i. This demonstrates, in slightly more formal terms, that the leading IR hamiltonian 
looks as (|1.1|) , but with a perturbation A^ cos (with dimension A^ = where $ is 

the dual of the free boson in the usual sense, and by dimensional analysis, A^ oc A~s. 



The same hamiltonian ( |1 . 1|) appears also in the problem of tunneling between edge 
states in the fractional quantum Hall effect 0. In that case, while it is Laughlin quasi 
particles of charge g = v (the filling fraction) that tunnel in the UV, the duality argument 
demonstrates that it is electrons of charge unity that tunnel in the IR. 

The duality just discussed is very useful qualitatively. It has however been used in the 
literature as a much stronger statement: namely that physical properties should exhibit 
an exact duality between the UV and IR fixed points under replacement of A by A^ and 
g by Why this should be the case was not explained, and it must be stressed that 
this is a highly non trivial result: the approach to the IR (strong coupling) fixed point 
is, in general, determined by a very specific combination of irrelevant operators coming 
with amplitudes that are all powers of Ad, so they all contribute equally significantly: for 



1 Space dimensionality does not play a crucial role here 
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instance, one expects that, in addition to the term Ad cos y terms X d Ad " 1 cos n^J 
should also appear (where = -), corresponding physically to multi-instantons processes, 
or tunneling of several electrons. Such terms might also be required as counterterms to 
cure the very strong short distance divergences of the IR perturbation theory. Clearly, the 
existence of these terms will destroy any hope of observing an exact duality, and one should 
not expect the duality argument to tell us more than the leading irrelevant operator, in a 
general situation. 

Nevertheless, an analytical computation of the mobility (the current) at T = and 
with an external force (bias) || has exhibited an exact duality between the UV and IR for 
the model adding up confusion to the whole issue. Something very special must be 

happening in that case - and indeed, the model is integrable. 

One of the purposes of this paper is to discuss why integrability gives rise to an exact 
duality for some physical properties - and also, to explain why this duality should not be 
expected for other properties. In discussing these questions, we will actually consider IR 
perturbation theory, and show how it can be made meaningful, again thanks to integra- 
bility. This has applications beyond the tunneling problem: as an example, we discuss in 
details the case of the resistivity in the Kondo problem. 

In the second section of this paper, we use the simple example of the Ising model with 
a boundary magnetic field to discuss how the integrable structure of quantum impurity 
problems gives a quick access to the full hamiltonian near the IR fixed point, which is 
essentially encoded in the reflection matrix, or the boundary free energy. We discuss the 
issue of regularization for the IR perturbation theory that arises in integrable models, and 
how one can use the renormalization group backwards in some cases. 

In the third section of this paper, we discuss how the IR action can be determined for 
the spin 1/2 Kondo problem and for the tunneling problem. We also sketch the result for 
the higher spin Kondo problem. 

In the fourth section of this paper, we discuss, as an application of IR perturbation 
theory, the resistivity in the Kondo problem - and determine, beyond the well known T 2 
order, its complete low temperature behaviour. 

In the fifth section of this paper, we finally discuss duality issues. We show that the 
anisotropic higher spin Kondo model exhibits a partial duality, manifest for instance in 
the following relation 





that holds up to analytical terms (odd powers) in H/Tb- We also show that the current 
in the tunneling problem obeys an exact duality 



These duality properties follow from the structure of the IR hamiltonians that is 
strongly constrained by integrability: within our "analytic" regularization scheme, they 
are made up of an infinite series of local (conserved) quantities (polynomials in derivatives 
of $), plus at most one non local term, which is A^cos ^/^$(0) for the BSG case, and 

\dS-e^~^^^+cc for the spinj Kondo case (where here S± are spin j — 1/2 operators) (all 
this within a well defined regularization scheme). As a result, thermodynamic quantities 
will in general exhibit a partial duality; the UV expansion in even powers of A will match 
the part of the IR expansion that is in even powers of A^, although there will also be other 
terms in this IR expansion due to the local conserved quantities. Some other properties 
turn out to be blind to the local conserved quantities however, and as a result exhibit an 
exact duality, like the DC current in the tunneling problem. 

In the first appendix, we determine the normalization of conserved quantities in the 
sine- Gordon theory. The second appendix contains some remarks about the Keldysh for- 
malism and analytic continuation. 

Some of the results presented here have appeared in short form in The methods 
we develop are related, although independent and different, to the series of works by 
Bazhanov, Lukyanov and Zamolodchikov ||, and also to the work of Lukyanov ||. Duality 
in quantum impurity problems has also been investigated by Fendley [[7], and by Fendley 
and one of us [§, from a more formal perspective. 

2. Getting the IR hamiltonian in integrable boundary field theories: the case 
of the Ising model. 

2.1. Some generalities 

We consider the Ising model defined on the half space x G [— oo, 0], y G [—00, 00]. We 
initially use a crossed channel or open string description, where euclidian time runs in the 




(1.3) 
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Fig. 1: Geometry of the problem. 



y direction, and we introduce the complex coordinate w = —y + ix. We restrict to the case 
of a theory that is massless in the bulk, and add up a boundary magnetic field h (Fig. 1). 

When h = 0, the fermions have free boundary conditions tpL = i^R on the boundary; 
this is the UV fixed point. When h — > oo, the Ising spins become fixed, corresponding to 
ipL = —ipR'i this is the IR fixed point. The question we wish to study is how the IR fixed 
point is approached. More precisely maybe, we want to be able to describe the Ising model 
at large values of h with a hamiltonian 

H = H m + 6H(0) 

where 7i expands in some powers of the inverse coupling constant 1/h. It is possible to 
gain some quick insight on what 5H should look like. - it is actually an expansion in 
odd powers of 1/h 2 . This is because the operator content of the Ising model with fixed 
boundary conditions can easily be extracted from conformal invariance considerations: 
with fixed boundary conditions on a cylinder of length L and circumference 1/T, the Ising 
partition function is simply the identity character (setting q = e _7r / LT ) 



1 

Xo=2 



oo oo 

n (i - q n+i/2 ) +u{ i+ i n+i/2 ) 



From this, it follows that the only available operators are of the form d p ipRd q ipR-\-{R — > L) 
(here, d stands for d w ). Up to total derivatives which do not affect the physical properties 
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of interest, we can restrict to i(>Rd n i(jR + (R — ► L), with n odd. Introducing the operator 
(the normalization is chosen for later convenience) 

0° 2k+2 = (-l) k+1 \ (: : + : ^df^L :) , (2.1) 

we thus expect 

oo ^ 

(0), (2.2) 

where the coefficients a 2 k+i have to be determined. Notice that in practice, the manipu- 
lation of expressions like ( |2.2j ) will give rise to extremely strong short distance divergences 
- the numerical values of a 2 k+i will only have a well defined meaning within a specific 
regularization scheme. 

For a general problem, such a computation would appear untractable. What makes 
it feasible in the cases we are going to consider in this paper is integrability. To see how 
this comes about, and pave the way for generalizations, let us describe the Ising model 
using massless scattering. In this simple case, we have massless R and L moving fermionic 
particles with energy and momentum parametrized as e = ±p = e^ 3 , (3 the rapidity. The 
mode expansion of the fermion operators is 



1>rW = J f^ /2 [uew(ePw)Z a (J3)+uew(-e?w)Z* R {l3)] 
M™) = J |^/ 2 [ueMe P w)Z L {(5) +ueM-e fi w)Zl (/?)], 



(2.3) 



where the Z are creation and annihilation operators obeying the usual anticommuting 
relations, uj = e l7r / 4 . The theory is defined on the half space x G [— oo, 0] only; as a result, 
the L and R modes are not independent. Because the boundary interaction is integrable, 
the fermions, in the crossed channel picture scatter off the boundary one by one, with no 
particle production, and one has Z R ((3) = R{(3)Z^{(3), R the reflection matrix ||. 

We will also use the direct channel or closed string picture, where euclidian time runs 
in the x direction. The mode expansion of the fermions is identical to ( |2.3|) , with w replaced 
by the variable z = w/i = x + iy. The Hilbert space is then the usual one for fermions 
defined on the whole line, and there is no relation between L and R modes. Rather, in the 
direct channel picture, the effect of the boundary is taken into account by the existence of 
a boundary state, which reads i 

ib>«E^ / Xl^K^- (3 B )zm)zmm, (2.4) 



2 We do not discuss the problem of the overall normalization of the boundary state in this 
multiparticle description - it is enough to recall that it is independent of (3b- 
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with K(J3) = R (f - (3). In the simple case of a boundary magnetic field considered so far, 
K((3 — (3b) = itanh ^ B 2 ~^ . The parameter is in general related with a typical energy 
scale associated with the boundary interaction, Tg = e^ B . In the case of a boundary 
magnetic field, Tg oc h 2 . In the closed string channel, we introduce the equivalent of ( |2TI| ) 



(2.5) 



2.2. The complete IR action 

Let us now discuss how the IR action can be simply extracted from the knowledge of 
the reflection matrix, or, equivalently, of the boundary state. 

To do so, let us keep working in the closed string channel, and consider the expression 
for the boundary state \B > further. The IR boundary state (fixed boundary conditions) 
is obtained as (3 b — » oo where K = i: 



n=0 



i s ^> « E h. I n ^ z lW)zm)\q > 



cm 

2n 



One can thus write \B) = B\Bm), where the operator B is defined in the multiparticle 
basis by 



bY[z*M)zm)\o>=U 



K{pi-(3 B ) 



Z2(A)^(A)|0> 



Let us expand 



In 



K{(3 - (3 B ) 



E 

fc=0 



2k + 1 



,(2fc+l)(/3-/3 B ) 



(2.6) 



Introduce then the set of commuting' operators X2/C+1 acting on the multiparticle states, 
with 



^2fc+l|/3l • • -Pn >C 1 ,...,C n = g 



(2.7) 



where C = L, R designates the chirality. One can then write 



|.B >= exp 



,fc=0 



2fc + 1 



e -(2k+l)f3 BZ 



2k+l 



\B m ). 



[2.1 



Of course, the I2/H-1 can be expressed in terms of the creation/annihilation operators, 
J 2fc+1 = J fg e (2fc+l)/3 [ZI((3)Z L ([3) + Z* R ([3)Z R ([3)}. Using the mode expansion of the 
fermions, one checks this coincides with dy(D2k+2, where 02fc+2 is defined in ( |2.5|) . 
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We can now write a reasonable conjecture for the hamiltonian ( |2.2| ) in the crossed 
channel - the reason why it is a conjecture only is because the exponential in ( |2.8|) is 
determined by the action on one particular state only, \Bm), and not in true generality (one 
can determine the action of the exponential on other states with "momentum" actually, 
but still, not on all possible states of the theory). Observe now that if Ti = TCm + 57i, 
the boundary state will generally read \B) = V exp — f_ dySH \Bm), where V is the 
(y) path ordered exponential. Using that the T 2 k+i form a set of commuting quantities, 
together with the fact that the 2 k+2 are self and mutually local operators, we obtain 
therefore 

oo _ 

H = H IR + J2 2^ e_(2 * +1)/3B0 2 fc+ 2(O), (2.9) 

k=0 

again up to total derivatives. 

Because the T 2 k+i form a set of commuting quantities, the perturbation of the IR 
hamiltonian in (|2.9|) is, formally, integrable. This is an expected result, since after all the 
flow from the UV to the IR fixed point is integrable, a feature that should be observed 
from both extremities - and provides an immediate check of ( |2.9| ). 



2. 3. The boundary free energy 

We now discuss the relation between the IR hamiltonian, and the boundary free 
energy. 

Consider thus the theory defined for x G [— oo, 0] and y G [0, 1/T], with periodic 
boundary conditions in the y direction. In the closed string point of view, the theory is 
thus defined on a circle instead of the infinite line, while in the open string point of view, 
it is now at finite temperature T. 

To compute the free energy, it is convenient to "unfold" the problem, so now x G 
[ — oo, oo], and the boundary interaction becomes an "impurity interaction" acting only on 
the R movers. Notice that a different unfolding is appropriate to study the vicinity of 
the UV and the IR fixed point; in one case, one extends the theory to x > by setting 
ijjR(x,y) — i/jL,(—x,y), while in the other one sets of course %j)n{x,y) = —^L^—x^y). In 
what follows, we discuss mostly the vicinity of the IR fixed point, and thus use the second 
folding. From the resulting "impurity" point of view, the IR fixed point is then just 
described by free R moving fermions. 



Introducing then, in the closed string channel 

1 JT 1 jT 

l2k+i = [ 2k+2 (z)dy = \ f dy: ^ R df +1 ^ R :, (2.10) 
Jo 1 Jo 

it follows from the expression of the boundary state that i 

oo 

/ = -T\ng IR - ^— -e-( 2fc+1 ^ 1/T (0|J 2fc+1 |0) 1/T , (2.11) 

k=0 

where gi R = i/t(0|-B/_r)i/t is the boundary degeneracy of the IR boundary state (actually 
independent of T |I(J) and |0)i/x denotes the ground state of the theory on a circle of 
circumference 1/T. 

The ground state on a circle corresponds to fermions with antiperiodic boundary 
conditions. Using the mode expansion of the fermions, and the expression of I 2 k+i as the 
sum of the 2k + 1 th powers of the energy, it follows thati 



1 oo 

1/T (0|J 2fc+1 |0) 1/T =-(27rT) 2/£+1 (0| £ (i + l/2) 2fc+ V-,-iM-+i/2|0) 

7 = — oo , . 

(2-12) 

= - 1 -(2irT) 2k+1 J2U + m 2k+1 - 

3=0 

The sum can be evaluated by ^-function regularization leading to 

1/T (0\l 2k+1 \0) 1/T = i (27rT) 2h+1 (l - C(-2k - 1), (2.13) 

(the same computation would give a vanishing result for even powers due to £( — 2k) = 0). 



For k = one gets — with c = 1/2; this is because X\ is nothing but the zero mode 
of the stress energy tensor § on a circle, Xi = {2nT) (Lq — By plugging the results 
( |2.13|) back in (|2.11j ), one obtains an explicit expression for /. 

Of course, the computation can be done in the open string channel as well; the ex- 
pression ( ^.9| ) for the hamiltonian leads to 

oo 2 

/ = -T\ng IR + £ ^—e-V k+1 ^(0° 2k+2 ) T . (2.14) 

k=0 



Q 

Here we have subtracted all extensive non universal terms. 

4 We use here the well known fact that the "conformal normal ordering" is related to the 
"operator normal ordering" by zeta regularization of the divergent parts. 

5 For the Ising model considered here, T zz = u : ^rOi^r : 
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Here, (0% k+2 )T follows from expression (|2.5|) evaluated in the multiparticle basis of the 
open channel: 

(0° 2k+2 ) T = {~l) k+1 \ J ^ 2k+1 ^(Z* R (f3)ZR({3))T, (2.15) 

and the notation (.) y designates the thermal average in the theory at temperature T. 
A standard thermodynamic analysis gives 

(0° 2k+2 } T = (-l) fc+1 T / d(3e^ k+1 ^ p((3) 



L 



oc 
oo 



i r°° i 

= (-l)^ r J_J^^ (m+m) T — m (2.16) 

1 ! J-oo^ d/31 + exp(e/T)' 

where, for free fermions, e = e 13 , and thus 2tt(p + p) = Ljjfe, a result that generalizes to 
interacting theories. By expanding the filling fraction, one obtains 

(0° 2k+2 )r = (~l) k+1 ^±ll T ^ + 2 ^ _ _L_^ C{2k + 2) . (2.17) 



To compare ( |2.13| ) and ( |2.17| ) recall the identities |TI 



((2k + 2) = V / ■ (-lfB 2k+2 - ((-2k - 1 



2fc+2 



2(2fc + 2)! v 7 " v " r " ^ 7 2A; + 2' 

where B n are Bernouilli numbers. Hence, as of course should be, ^ 2k+2 )T 
T 1 / T (0|X 2 fc+2|0)i/T5 an d we find the same expression for the impurity free energy |12|| . 



Using the thermodynamic expression for the integrals of motion (the first equation of 
( [2.160 ), we obtain an alternate formula for the impurity free energy: 

2 J2k+l)(B-B B )r i\k+ide 1 



f = -T\ng IR + [ ph+W-p^yk+i 
yiH J 2?r 2k + 1 y ' 



d01 + e<M T 

(2.18) 



-Tln 9uv -T ! f- 1 Ul + e-W) 

J 2tx cosh(/5 — pb) ^ ' 



where we used the fact that / ~ —Tin guv (resp. / ~ —Tlngm) as (3b — » — oo (resp. 

~~ ^ oo). This last expression coincides with a well known formula obtained using the 
thermodynamic Bethe ansatz. It reads as well 

f = -Tlng uv -T f ^--^ In R(/3 - B ) In (l + e~^ T ) , (2.19) 

a result that follows directly from the form of the boundary state, and the manipulations 
in flTT6D . 

In the foregoing paragraphs, we have thus showed how the IR action could be extracted 
from the R matrix, and how it was closely related with the boundary free energy. 



9 



2.4- Flowing "back" from the IR fixed point. 

The previous analysis shows very clearly how the TBA results are directly connected 
with an IR description of the flow; in fact, the free energy provides an immediate reading of 
the complete IR action (that the impurity free energy has to do with conserved quantities 
was observed in the earlier papers on the subject already, see [fl3]). It is important to 
realize that all this works for a particular regularization scheme, involving dimensional 
regularization and (or) contour deformation. This is somewhat obvious since the integrable 
approach does not involve any length scale that could act as a cut-off. The quickest way 
to see this more explicitely is to consider for instance the quantity X\. Using the mode 
expansion, \0)i/t is clearly an eigenstate of Xi, and thus 

1/T (0|(X 1 ) P |0) 1/T = ( 1/T (0|Xi|0) 1/T ) P = (-^) P - (2-20) 



To write (|2.20| ), we have used an operator formalism, which, in fact requires "time ordering" 



- here ordering along x. In other words, in ( 2.20 ), the divergences have been regulated by 
slightly displacing the p contours of integration. 

The effect of this displacement can be seen by using the fermion propagators and 
Wick's theorem. Aside from the term involving the (non vanishing) average of X\ on the 
circle, contractions contribute integrals with strong short distance divergences, the simplest 
one being 

l/T 1 

dyidy 2 - — — -4. 

[sin7rT(yi - y 2 )\ 

If one evaluates this integral by displacing the contours and using the residue theorem 
(together with the periodicity of the integrand), one finds indeed a vanishing result, because 
the integrand has a vanishing residue at the origin. 

Equivalently, in dimensional regularization, one considers the more general integral 
where the power is a number a (a = 4 here) and one computes the integral in the domain 
of ck's where it is defined. This gives 



IT 1 



TV 



2 1+a r(2 - a) 



(l-a)I*(l-f) 

One then continues analytically to a = 4 - and the last expression vanishes again, this time 
due to the double pole (in a) in the denominator. This generalizes to all the other integrals, 
so that the dimensionally regularized value of X\ is contributed only by its average, ie the 
result (EE2CD . 
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In general, integrals involving local conserved quantities (all the ones in the Ising model 
are of that type) will be regulated by operator methods or contour displacement, since 
there is no readily available parameter to perform continuations (the prescription is not 
ambiguous thanks to the commutativity of the conserved quantities). Integrals involving 
non local conserved quantities will be regulated by continuation in the parameter g. We 
refer to this scheme as an "analytic" regularization. 

Using the previous ideas, it is clear that we can solve the problem of the most general 
perturbation of the IR fixed point. For an action 

oo 

n = Hi R + Y t hk+iOS k+2 , (2.21) 

fc=0 

the boundary free energy simply reads 

/=-Tln SlH + / g g Wl e<^(- 1 )^! T _L^, (2.22) 

the integrals themshelves being evaluated in ( [2.16| ),( pTT7| ). 

As the temperature is lowered, ie when one considers this system at larger and larger 
scales, one simply flows to the IR fixed point, as physically expected, since all the operators 
^2fc+2 are irrelevant near this fixed point. As the temperature is increased, ie when one 
considers the system at smaller and smaller scales, or tries to "flow back" , what happens 
generically is that no fixed point is reached; rather, the amplitude of all the terms becomes 
bigger and bigger, as expected for irrelevant perturbations. The cases where one flows back 
to an interesting fixed point are the ones for which the series in ( |2.22|) defines a function of 
T which, continued beyond the radius of convergence, has a finite T — > oo limit. Though we 
do not know any definite mathematical statement about that question, it seems clear that 
these cases are extremely rare. For instance, the choice &2fc+i = 2fc 2 ) _ 1 e~^ 2fc+1 - ) ^ B guarantes 
a flow back to the free fixed point, but any perturbation that differs, even infinitesimally, 
from this one by a finite number of terms, will not flow back to the free fixed point at all. 

A quick way to build an IR hamiltonian that has aT^oo limit is to multiply the 
reflection matrix by a CDD factor. By flowing backwards, one finds in this case that 
the difference guv ~ 9ir is increased by a term In \f2, corresponding presumably to the 
appearance of additional boundary degrees of freedom in the UV. 
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3. Approach to the IR fixed point for the spin 1/2 Kondo model and the 
boundary sine-Gordon model 

3.1. The Kondo model 

The previous structure generalizes in a slightly more complicated form to the case of 
the spin 1/2 Kondo model with action 

1 t>0 poo 

U= -\ dx dy [(d x $) 2 + (II) 2 ] +X S_ e V^(0) + >s+e -V2^(0) 5 (31) 

where S is a spin one half operator (A is assumed positive in what follows). The boundary 
interaction is integrable, and the same manipulations we carried out for the Ising model 
can be acomplished here too. 

Instead of describing the bulk with massless fermions, we use massless L and R moving 
solitons and antisolitons. Parametrizing their energy by a rapidity e = ±p = e@, these 
particles have factorized scattering, the LL and RR scattering being given by an S matrix 
which, as a function of the rapidities, is the same as the S matrix of the bulk sine-Gordon 
model, Sll = Srr = S SG while the LR scattering is trivial. The solitons and antisolitons 
scatter off the boundary one by one with no particle production, and the R matrix is given 
by 



4 J (3.2) 



R% = 0. 



In the so called repulsive regime, that is for g > |, there are no bound states, and the 
soliton and antisoliton are the only particles in the spectrum. The boundary state can be 
written in a form similar to (^.4j) 



fo J ^<-<^ i 271 e^± (3.3) 
xZ£ i£i (/? 1 )...Z£ i£ J/? n )|0), 

with, as for the Ising case, K{j3 — 0b) = itanh ^® ■ As a result, introducing (C denotes 
the chirality, C = L, R) 



WIA-A>c ll£1 ... = ^ ^ e (^i)/^ |A.../3 n ) Cii£i . 



(3.4) 
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we can write 



| B) = exp 



E 

.k=0 



-2 



(2fc + l) 



■(2k+l)0 B J 2fc+1 



A 



2fc+l 



|Sj fl >. 



(3.5) 



The coefficients A2/C+1 in (|3.4| ) will be adjusted for later convenience. 

Indeed, a new difficulty arises here when one wishes to reexpress the set of commuting 
quantities T2k+i in terms of local fields. As far as we know, this question was first addressed 



quantitatively in ||12| , where the first few conserved quantities were studied numerically 
using the TBA. The following analytical expression was obtained in unpublished works by 
the present authors, as well as by Al. Zamolodchikov ||14j| , and probably by a few others 
too. A derivation is presented in the appendix for completeness; to our knowledge, it has 
never appeared elsewhere, though the technique is hardly original. 

To proceed, we need to chose some normalizations. We first introduce the twisted 
stress energy tensor 



2tt 



-27t:(dcj>y :+i(l-g) x l— d 2 




(3.6) 



where (p = <J)r is the right moving component of the boson. The central charge correspond- 
ing to this tensor is 



1-6 



(1-9)' 
9 



(3.7) 



A set of commuting quantities is then obtained by integrating successive powers of this 
stress energy tensor. We define (the 2tt normalization makes subsequent formulas simpler) 

1 



o 2 = -r (t zz + Tzz) 



o 4 



±-(;T* z + R^L) 



1 



c + 2 
12 



: T zz d T zz : +R 



(3.8) 



The normalization is such that O2/C+2 goes as |(— l) fc+1 (27r) fc (9^) 2fc+2 . We then define 
%2k+i = /_„ 02k+2dy. Let us stress that these quantities commute at the conformal point 
only (in fact, of course, their left and right components independently commute). In the 
massive sine-Gordon model (ie with the bulk perturbation cos2y / 2~7r^$ in our notations), 
there exists non chiral deformations of these quantities that still commute, and act as 



sums of odd powers of momenta on the (massive) multiparticle states [ 15| . In the massless 
scattering description we are using here, one considers the free boson as the limit of the 
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massive sine-Gordon model, and the particular chiral quantities T 2 k+i are singled out, 
which act again as sums of odd powers of momenta on the multiparticle states. Of course, 
there are more conserved quantities right at the confer mal point, but they do not seem 



to have any simple meaning in terms of rapidities - see next section however, and [T^] for 
more details. 

With this choice, one has (see the appendix) 



A 



2k+l 



(fc+1)! 



"(2A;+l)g" 




1 


) 


[ 2(1-9) _ 




L 2(1-9) J 





2fc+l 



9 




r 


"(2/c+l)" 


|.2(l-fl) J 


^ 2k+l 




[2(1 — ff) J 



(3.9) 



In the following, we will also need the relation between the parameter (3b of the R matrix 
and the coupling A in the action of the Kondo model. This was determined in ||, and 
reads 

T B = yf- 9) \ [AT(1 - g)] l/{l - 9) . (3.10) 



2(1-9) 



From ( |3.5| ), we then obtain 



2fc+2' 



(3.11) 



with 



J2k+1 



2k + 1 



-(2k+l)f3 1 



k=0 



l + 2fc 
A !-3 



A 



2fc+l 



n k+i r 

y 


'(k + 1/2)^ 




(l-g)(k + l)l T 


1 + (fc + 1/2)^2 


9 



l+2fc 

[r(i-^)]-— 



(3.12) 



and O2/C+2 follows from the expression for 2 k+2 by replacing z by w, and multiplying by 
an overall factor (— l) fc+1 : 



o° 2 
o° 4 



4:7V 



4tv 



1 ' . 7-3 

47T l " WW 



c + 2 
12 



(3.13) 
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The foregoing results essentially coincide with those in [H . Our route is quite different 
however; in particular, the form of the R matrix or the normalization of the integrals of 
motion are not used at all in H , where, instead, a functional relation approach is developed. 

The Kondo model is a very interesting physical example from the point of view of 
the IR perturbation theory. For any value of g (which physically corresponds to the 
anisotropy), the IR fixed point is always the same (see [[16] for details and references). To 
get back to a g dependent UV fixed point, one needs to perturb the IR fixed point by the 
same family of operators (stress tensor and the like) but with coefficients that depend on 
g: it is only through this fine tuning of the coefficients that different flows can be obtained. 
For a given g, the free energy for an arbitrary IR perturbation - that expands on the 
conserved quantities - has an expression similar to what we wrote in the Ising model. 

The foregoing analysis could be generalized to the regime where the associated bulk 
sine-Gordon model has bound states, that is g < \. The final expressions involving 
quantum fields, for instance ( |3.11] ), would not change; they are expected to be analytical 
in g, a result that can easily be checked using the method we explain below. On the 
other hand, expressions involving scattering quantities would change . Here, we would 
like to make a remark concerning ( |3.4|) . Because in the scattering the numbers of solitons 
and breathers are independently conserved, one expects in general a result of the form, 
introducing the color for particles (e = 1, . . . , m, . . . for breathers, e = ±1 for solitons 
antisolitons) 

l 2k+1 \p 1 ,...,p n ) Cl , ei ,... = \ (^£\2k+i, ai e {2k+1) ^ |/3i,...,/3 n ) Cl , £l ,..., (3.14) 

where the A2fc+i,a are a priori all different. The determination of these factors is an 
interesting problem by itself. It can be quickly solved if one observes that the formula for 
the boundary state ( |3.3| ) immediately generalizes to the case where breathers are present 
in the spectrum, by using the m-breather reflection matrix 

tanh(^-f^) 

Rm = 7 -4r. (3.15) 

tanh(^ + f^) 

Expanding i^lni? in (odd) powers of e 13 , one finds 



-— ln J R m = 4y e -( 2fc+1 ^sin 
1 <® fro 



(2fc + l)g 
run — ; — 

2 1-0) 



ld_ 

idj3 



oo 



(3.16) 



l nj R = 2^(-l) fc e-( 2fc+1 ^. 
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By putting these expansions in the formula for the boundary state, it follows that the ratio 
of normalizations of conserved quantities is the same as the ratio of the odd powers of e 13 
in ( p.!6| ), that is 

(2k + l)g~ 



^±I^ = 2(-l) fc sin 

A2fe+1,± 



mm- 



(3.17) 



2(1-3) 

The rest of the arguments follows with minor modifications. 

The key feature of the spin 1/2 Kondo problem is that the IR fixed point is approached 
along the conserved quantities C*2fc+2 °f even dimensions. The situation is more interesting 
for the higher spin case, or the boundary sine-Gordon case. 

3.2. The boundary sine-Gordon problem 

The previous structure generalizes in a slightly more complicated form to the boundary 
sine-gordon model 

i M poo 

H = - dx dy [(d x $) 2 + (IT) 2 ] + 2A cos v^^$(0). (3.18) 

J — oo J — oo 

The boundary interaction is integrable, and the same manipulations we carried out for the 
Ising model can be acomplished here too. While technically more involved, the general 
spirit is very similar, so we will restrict ourselves to the salient features. 

The quickest way to proceed is to restrict to the attractive case of the associated bulk 



sine-Gordon model, g = 1/integer, and to consider the reflection matrices in that case ||17|| : 

lAinR -off 1 \k p -(2k+i)p sin [mn(2k + l)g/2(l - g)} 
id(3 mKm ^ } sin [7r(2fc + 1)0/2(1 -g)) 

-4^HRt ± R -) = y(-l) k e- 2k ^-^ tan ibr^— ^ (3.19) 

id/3 g k~i 9 

1 00 00 — (2fc+l)/3 

[ ~ 9 J2(-l) k+1 e- {2k+1)l3il ~ 9)/9 + ' 



9 fo ^ Q sm[n(2k + l)g/2(l-g)Y 

The boundary scattering is non diagonal in the soliton antisoliton basis, but it is diagonal 
for the symmetric and antisymmetric combinations, which scatter with the amplitudes 
± Rt. = R±. The bulk scattering is diagonal in either basis. 

We also recall that the relation between the coupling constant A in the action and the 
rapidity (3b is modified in the case of the boundary sine-Gordon model, reading then 



T B = (2smirg) 1/(1 - 9) ^ g) \ [AT(1 - g)} 1/(1 ~ 9) . (3.20) 



2(1-3) 
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By using the reflection matrices, and following the same logic as before, we obtain 
immediately the coefficients of all conserved quantities for the hamiltonian near the IR 



fixed point. We have now an expansion similar to (|3.11| ), but with the coefficents 62/c+i 
replaced by 

(-l) fc 1 1 

C2/C+1 = — ~ : — j—7771 — r~T\ — ToTi vi 2 fc+i b 2 k+i, (3-21) 

2 sin [n(2k + l)g/2(l-g)] ( 2s in7r#)^) 

where the prefactor is just the ratio of the coefficients of odd powers of e-* 3 in ( ETIgD and 
( |3.16| ) (of course, this ratio is the same for the breathers and the soliton antisoliton R 
matrices) , plus an additional power of 2 sin rrg arising from the difference between ( |3.2Q ) 
and (UUCP . 



It is well known |T^] that, at the conformal point, the (chiral part) of the quan- 
tities l2k+i commute not only together and with the integral of the perturbation 

dye^ 1 ^ 8 ^ 9 ^ , but they also commute with the "dual" of the perturbation e ±l ^~^^ . 
When the perturbation is turned on, a deformation of these quantities turns out to still be 
conserved, guaranteeing the integrability of the flow. This conservation is true all the way 
to the IR fixed point, where again the purely chiral quantities are conserved, by conformal 
invariance. It follows that, if one investigates the conservation perturbatively near the 
IR fixed point within the dimensionnally regularized scheme, the only operator that can 

2tt, 



be added to Hir, besides the O2/C+2, is cosy-^$. Here, $ is the dual of the field <f>, 
3* = 4>R — 4>l, and we used that <fin = —<Pl at the IR fixed point. By dimensional analysis, 
its amplitude goes as A _1 / ff . The exact amplitude follows from eqn (6.18) and (6.20) in 
0. One can thus finally write 

[2~K °° l+2fe 

H = H IR + 2\ d cosJ—$ + J2 c 2k+i*~ 1 ^OZ k+2 , (3.22) 

' 3 k=0 



where 



A d = — r 

2ng 



(;) 


[<7r(<7)l 




2tt 



9 X~9. (3.23) 



Observe that the R matrix elements for breathers expand only on odd powers of - 

eigenvalue on the breather states, a result of their charge neutrality. 
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3.3. The Kondo model with higher spin 



It is necessary to still generalize the previous arguments slightly, to take into account 
the Kondo model with higher spin. The structure is actually very similar to the spin 1/2 
Kondo and the BSG case. 



The UV Kondo hamiltonian reads as (|3.18|) with now the spin in a spin j representation 
of U q sl(2), q = e 1 ^ 9 . One finds that the hamiltonian near the IR fixed point reads 



T~t = Hir + Xd 



S-e' 



— 5(0) . a ~ 

9 w + S+e 



l-(O) 



(3.24) 



k=0 



where s is in the representation j — |, Xd oc A and the coefficients d 2 k+i could be 
determined using the same method as before (see section 5 for more details) . 



4. An application: the resistivity in the Kondo model. 

A good testing ground for the previous considerations is the isotropic Kondo model, 
where the strong coupling behaviour can be probed by experiments at low temperatures. 
The most interesting quantity in that case is of course the resistivity, for which no closed 
form results were available so far, besides the T 2 term that follows from Fermi liquid theory 
fl9fl (attempts to compute p with the Bethe ansatz have failed, partly because it is truly 
a three dimensional quantity). The method we have developed in this paper allows us 
to make an important progress on that question: short of getting p in closed form, we 
can at least compute it perturbatively near the strong coupling fixed point, now that we 
know the exact structure of the hamiltonian. This allows us to go beyond the Fermi liquid 
approximation, and evaluate p as a power series in T 2 at low temperatures. 

In order to study the resistivity, we first need to go back to the 3d formulation of the 
system with electronic anihilation operator 

* (f) = / 1075°***®' (zL1) 

where we suppressed the spin indice for simplicity. As usual, since the Kondo interaction, 
when the impurities are dilute, is assumed to be with only one impurity, we can consider 
only the s-wave component of that operator around the Fermi points 

tt(f) = _^_[e**'^fl(r) - e-^ r Mr)} (4.2) 
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with r > and we used right and left one dimensional moving field. This decomposition 
implies V'-l(O) = i/jr(Q). In the interacting theory, only the s-wave parts of the three 
dimensional Green's function will be affected, moreover, only the LR and RL components 
of the dimensionally reduced model are affected by the interaction. This leads to the 
following form for the three dimensional interacting Green's function (for the spin up field 
for example) 



G(u M ,n, r 2 ) - G (u M , n - f 2 ) 
-1 

g ?r 2 rir2 

+e ikF{ - ri+r *\G RL (u M , ri, r 2 ) - G RL (u M , r u r 2 )) 



e- ik ^+ r *\G LR {u M , n, r 2 ) - G° LR {u M , n, r 2 )) 



(4.3) 



with the superscript denoting the free Green function. As we will see, this is the quantity 
we need to compute the resistivity. The interacting LR (resp. RL) Green's functions are 
defined by 

G Cl c 2 ("M,r) = - dy e^WcAwWczM), (4-4) 

J-P/2 

with Ci indicating the chirality. As an example, we haveH 

,/3/2 



IUIMV 



dy- 



p/2 § sin |(-y + zr) 



(4.5) 



where we have used the fact that r > and that in the UV 

(<Mwi)</>l(w 2 )> = —s , — — • 

^m^wj -w 2 ) 



(4.6) 



In the IR, the only difference is the boundary condition which will result in a change of sign 
in the propagator. When we put everything back into the three dimensional expression for 
the Green function, we get (at the IR fixed point) 



G ir (ljm, n, r 2 ) - G°(u M , ft - f 2 ) 

= 1 e -!tF(ri|r 2 ) e WM(ri+r2)^_ UM j _ (n +r 2 ) e ~^M (ri +r 2 . 

2nrir 2 L 

= G (ujM,f 1 )T(iv M )G (ujM, -f 2 ). 



(4.7) 



The fermions operators have an extra \/2-k in their normalisation here. 
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Following the arguments of [20], for a dilute array of impurities of densities nj the lowest 
order correction to the complete Green function takes the form 

G(u) M ,ri, r 2 ) - G°(uj M ,r 1 - r 2 ) ~ rn \ d 3 ; fjG c \uj M , fi - n)T \uj m )G z \u m , n - r 2 ) (4.8) 



(4.9) 



Summing over multi-impurity terms, the self-energy takes the simple form 



where higher orders in m are neglected. The retarded self-energy is found by the analytical 
continuation iuiM — ► oj + it] leading to 



(4.10) 



v is the number of spin per channel, we have reestablished its dependance at the end 
since it only amounts to a factor of two (separate spins contribute the same). This is 
the expected result at the IR fixed point for the one channel Kondo model. Finally the 
resistivity follows from the Kubo formula for the conductivity 



1 



cr(T) 



dn 
de k 



(4.11) 



p(T) ~ v ' 3m 2 J (2tt) 3 

with the single particle lifetime defined by 1/r = —2Im'E R . The dispersion relation e k = 
vpk has been linearised in that limit. 

All this discussion was done using the fermions but to continue and understand how 
to get away from the IR fixed point, we need to use our earlier results. To make contact 
with our previous discussion of the Kondo model, we need to bosonise the system. This is 
done using the rules 

^Pl/rA^ V) « e^^V^CMO. (4.12) 

Notice that we have reestablished the spin dependence, \i =|, j, since this will be crucial 
in the following. At the UV fixed point we have ^£,^(0) = ^^(O) but since we are 
interested rather in perturbation around the IR fixed point, we impose the conditions 
V'-RjmW = ~i J L,Li(0) for the IR correlators. This leads to the RL (LR) bosonic propagator 



(<I>r,u[ w i)<I>lA w 2)) 



1 , P . 7T 

— In - sin — (wi - w 2 ) 

47T 7T p 



(4.13) 



which translates in the correct fermionic propagator when using the bosonisation rules 
given above. Although we are interested in computing the Green function of spin up 
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fields, for example, the integrable description is rather in terms of the spin and charge 
densities, ie introduce 

f (4.14) 

<Pc = + < f>i)- 

In terms of these fields, the interaction at the boundary only involve the spin field and is 
given by the hamiltonian written in the previous section. The charge field remains non- 
interacting. The perturbation around the IR fixed point is described by the hamiltonian 

oo 

U = H IR + b2k+i\-^0° 2k+2 (4.15) 

fc=0 

where all the couplings and operators have been given in section 3; the boson field in the 
latter section coincides with S here. On the other hand, if we look at the bosonisation 
of the Green function for the spin up field, we observe that there will be contributions for 
each field 

and when computing the interacting left-right Green function for example, the charge 
sector will be completely decoupled, ie 



(•••) — (•• ^charge X (' ' ') 



spin 



Only when doing the Fourier transform will the charge part contribute. Let us proceed to 
the computation to show this more explicitly. 

The isotropic case (g = 1) leads to some simplifications in the previous expressions, 
leading to the identification 

(l + 2fc) 1 

W^ Tt(fc + , )( , + 1) T B "+-'. (4.17) 

The coupling Tb now is identified with the usual Kondo temperature Tk (up to a nor- 
malization that is a matter of convention, and will be decided later) and the contribution 
of each operator is determined through these relations. This provides the information 
necessary to compute higher corrections to the resistivity from the IR fixed point. 

Up to order T^ 2 the contributions are exactly the same as the ones found previously 
since only one operator, the energy momentum tensor, appears to that order. It is at the 
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third order that the non-trivial approach to the fixed point will be needed since the second 
operator 4 will be involved. First let us proceed to reproduce results found before for 
the two first orders using our bosonised formulation. To first order, the leading irrelevant 
operator is (with the proper normalisation) 

1 1 (2ttT 2 ) 

- — (T ww + Tu,u,) = -[: (d w <p s ) 2 : + : (<9™</> s ) 2 :] 



An '2 24 

■ (3 



(4.18) 

_(2ttT 2 ) 
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where we have used the fact that the operator is inserted at r = to get the last line. 
The constant is a disconnected term that gets cancelled when dividing by the partition 
function to evaluate the correlator: we can thus forget about it in what follows. Inserting 
([4.1 8|) in the correlator (of the relevant RL or LR components) we get the lowest order 
contribution to the one dimensional propagators (again for the spin up field for example) 

— / dydy' e to wv/ e ±iv^F[^ /i ,c+0R/i i ,.](ri,») . (d ,0 s ) 2 : x 
T B J-p/2 ' ' (4.19) 

Xe TiV2n[<l>L/R,c+<l>L/R,sKr2,0)\ R 

with the subscript IR meaning that we evaluate the propagators with respect to the IR 
action. Note that the contribution from the charge boson decouples and the perturbation 
only affects the spin sector. Again let us write explicitely the RL component: we have for 
the first correction 

-1 r^ 2 ^-sm^(w 1 -w 2 ) 

5 W Grl = / d y dy ' e ^MV * £ii V. 

4nT B J_ p/2 [f sin f(w> - w 2 )} 2 [^ sin f («/ - (42Q) 

J-B 

with {(wm) the step function. This leads to a correction of the self-energy of the form 



to 

2 + i — 



(4.21) 



2tiv 

which is the expected form. The correction is real and does not contribute to the con- 
ductivity or the life time. To get bona-fide contributions, we need to go further in the IR 
perturbation theory. To next order, the conserved quantity Oi will contribute again but 
the higher quantity O4 will not yet give a contribution. So to second order, we have 

2 fP/ 2 piuMV 

-=2 / dydy'dy" 



T b J -0/2 [£ sin I ( Wl - W2)] 1 / 2 (4.22) 
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where we already contracted the charge part. Using the relation 

: OV0(O, y')f :: (0 y „0(O, y")f :=: (<V0(O, y')) 2 (<V0(O, y")) 2 

' ' :M0,i/))(V^/)): 

(4.23) 




^sin^y'-y")] 2 

we get three contributions to the second order, two of which are divergent. The regular- 
isation of divergences here is done by analyticity, as explained in section 3: we slightly 
modify the contours of the y" integral, and move it by id in the complex plane. The 
integrals are then done by simple residue evaluation. Usually there could be a dependance 
on the way the contour is deformed but this disapears here since the operators commute 
with each other (there is no simple pole in their OPE). The last term in the expansion 
has no frequency dependence and the explicit evaluation (using our prescription for the 
regularisation of the divergence) gives zero. The first contribution has the form 

S^G RL = - -2— / dydy'dy"e^y sin -{w, - «; 2 )] 3 x 

X< [fsmfK-^i)][fsmfK-«; 2 )]} X ^ 

1 



[£ sin |K - Wl )][% sin %(w" - w 2 ) 

which contains no divergences and can be evaluated straightforwardly by the method of 
residues. The integral over y',y" leads to 

. [cos %(wi - w 2 )} : 

d ye luJMy 4 — - 

P B J-p/2 [f sin|(«;i -w 2 )Y 

and evaluation of the integral gives 

in 

^2 
B 

The second contribution, which has divergences, takes the form 

-i r p/2 i 

S^G RL = \ dydy'dy"e lLJMV - R x 

^Tli_t3/2 \^%v'-y")? 

^sin|(iui - w 2 ) 



5^G RL = -==r / dye lU}My - n ^ — , (4.25) 



5 (2a) GflL = ^e(, M )e-- (ri+r2) [4 + (vrT) 2 ]. (4.26) 



(4.27) 



x 



[f sin - u) 2 )] [f sin - «^)] [f sin - w 2 )] [f sin § («;" - i^)] 
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and evaluating by residues leads to 



S (2b) G RL = ^e(u M )e-^+^(2uj 2 M - 2(ttT) 2 ). 



(4.28) 



So the total contribution to second order to the retarded green's function takes the form 
(once we analytically continue to real frequencies) 



2tvu 



00 1 

2 + i- 



- (3^ + W ) 



(4.29) 



and as expected we have a universal function of (uj/T k ,T/T k ). The previous results did 
not require any information about the other operators but at third order, the operator 



CO = — - [• T 2 • -\- • T 2 - •] 

*-^4 [• - 1 ww ■ ' ■ ww •! 



= n 



1 



: (d y <p s ) : -— : d y (j) 8 d y (j) s : 



(ttT) 2 fa x x2 , 3(ttT) 4 



(4.30) 



80n 



needs to be taken into accout: it comes with the coupling l/(37rT|) in the hamiltonian. 
Using the relation 



■■(dy 



: : e 



— i\/~2ir<j) L {w2) 



n—p 



X 



(4.31) 



we get, using the residue theorem, the contribution, 

S {3a) G RL = -^e(u M )e-^^+^[6(tu M ) 3 + 6tu M (nT) 2 ] 



(4.32) 



There is also a contribution from the leading irrelevant operator when expanded to third 
order, which reads 



^ b) G RL = J^ e (u, M )e- u » ^+^[3z^ M (7rT) 2 + 5(iiv M ) 3 ] 



(4.33) 



At this order, the contributions are all imaginary and we need to go the the next order to 
get non trivial contributions to the resistivity. At fourth order there are two contributions, 
one coming from the leading operator only, C?|, and another, from the combination of 
the leading and next to leading operators, O2O4. The computation are analogous to the 



24 



previous ones, but more tedious. The final result for the retarded self energy up to fourth 
order is 



u) = - 



2nu 



2 + €-^f( 3 - 2 + ( " T)2 )- 



5 3 



+ 



+ 



12 24tt 

35 7 
+ 



to 



B 



192 32tt 

11 3 

+ 



192 32tt 



ttTV 
Tb~) 



4 ' 8n J T B \T B J 

2 



7rT\ 



19 



96 16tt 



ttTV 



(4.34) 



Using the Kubo formula this leads to our main result for the resistivity (which we computed 
to sixth order) 

2 



P(T) 



3n; 



(7TVfv) 2 

( 47 



4 \T B J 



13 3 
+ 



+ 



: 



1 



53 



240 20tt 



TFT 



(4.35) 



10080 8tt 336tt 2 J \T B J 
In the following figure we compare this result with the numerical renormalisation group 



method [^TJ. The definition of T B is related to the usual Kondo temperature through a 
simple factor T B = —Tk- 
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The agreement is quite good considering that there is no fitting parameter. The Pade 
approximants were found to be very stable, and give a control of the curve p(T) all the 
way to T m Tk, which is right in the crossover region. It is thus clear that our method 
provides a good analytical understanding of the strong coupling resistivity. 



5. Another application: duality 

The general structure of the IR hamiltonians is given by a set of local conserved 
quantities, plus at most one non local conserved quantity. This implies some duality 
properties that we now discuss. 

5.1. Duality in Kondo with higher spin 



The main thing about expression (|3.24|) is that it contains only one type of exponential 



@. Qualitatively, this is a consequence of integrability: it is natural to expect the trajectory 
to appear integrable both from the UV and IR fixed point; on the other hand, theories with 
several harmonics are generally non integrable - therefore, only one harmonic can occur. 
Quantitatively, this leads to a very strong similarity of the physical properties expanded 
near the UV or near the IR fixed point, after replacement g — > in particular, quantities 
that are "blind" to the integer spin conserved quantities, if any, will exhibit a complete 
duality symmetry between the UV and the IR. 

To discuss the matter further, let us compute the boundary free energy at vanishing 
temperature and with an applied field 2HS Z (S z taking values j,j — in the 

representation of spin j). We introduce the quantity e2j defined by 



£ 2i 



du sinh 


[(*= 




TTLO 

2 


2yr sinh 


g 7TLJ 

1-9 2 





e»-2U(2j-l). (5.1) 



Here e is the quantity defined in eq. (6.9) of with the conventions of the appendix 
(M = 2, h = 1,V = 2U, e = 1); one has in particular, e\ = e of @. The parameter V 

7 Here, we should stress that this is a regularization dependent feature, that holds for our 
dimensionnally regularized approach. Ohter exponentials would appear as counter terms in other 
approaches. 
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is related with the field by H = gV. One can then establish, from the well known TBA 
formula in the limit T — > 0, that (this generalizes slightly |22^ . See also 0,|§) 
d(3 1 



2n cosh(/3 - (3 B ) 



e 2 M 



dw 



1 



(5.2) 



2tt 

In this formula, 



2 cosh 2g 



sinh 



9 7TOJ 

1-3 2 



cuia; — i) 



V(2j - 1). 



G_(uO 



2tt 




r[iw/2(i-^)] 



sr r[ W9 /2(i- 5 )]r[i/2 + !W /2] 



(5.3) 



and A = | In + 2(i- g ) In g. To compute /, we close the contour in the upper half plane 
when A > (3b- The only poles are those at u = 2(1 — g)ni, n a positive integer. The UV 
expansion of / follows 



-2(l- S )nA 



n=l 



sin2n7T(7 r (— ng) T [3/2 — n(l — g)] \e 



T, 



B 



2n{X-g) 



(5.4) 



2jVg. 



We now recall the correspondence between the cut-off A and the physical field in that 
case e A = V • Using that the field coupled to the impurity in the Kondo problem is 
H = gV ', together with the correspondence between the bare coupling A and Tg, one has 

e A He A 

(5.5) 



Tb [XTil-g)} 1 / 1 - 9 ' 
This allow us to rewrite the free energy in the form 



/ 



9 



E 



n=l 

- 2jH. 



nnl 



sm2mrg T (-ng)T [3/2 - n(l - g)] H^a 



AT(1 - g) 



2n 



(5.6) 



When A < (3b on the other hand, we close the contour in the lower half plane. There 
are now two types of poles: the ones at ui = — (2n + l)i give contributions to the free 



energy of the form ( ^ j 



V 



while those of the form ui 
(which we will refer to as "non-analytic" ) 

(_l)2ry sin ( 2 j - l)mr/g 

n=l 



-2ni—£- give the contribution 

9 to 



non— analytic 



nn\ 



Ar(i-g) 

H l ~9 



2n 
3 



sm2nn/g 
H 



n/g)T [3/2 -n(l-l/g)} 



(5.7) 



(2j - 1) 
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From this we deduce the relation (to be used in (|3.24| )) 

9 \~9, (5.8) 



sin - / i 

A rf = — - 



Kg \gj L n 

together!! with 

f(j,\,H ig ) = f(j-± A d ,-,-Y (5.9) 
V 2 g gj 

where the equality holds up to analytical terms (odd powers) in H/Tb- 

This duality has an obvious physical origin. We can compute the free energy near the 
UV fixed point perturbatively in powers of A, or near the IR fixed point perturbatively in 
powers of A^ oc \~ l / g and in powers of A -1 ^ 1-5 -*. The firs type of terms comes from the 
Kondo type interaction near the IR fixed point, that looks formally like the one near the 
UV fixed point, but with the replacements j — > j — 1/2, g — > 1/g and H — > H/g. It is 
interesting to discuss the later replacement in more details - the physical interaction near 
the UV and IR fixed points of course does not change, it is always 2HS Z . However, to 
take this into account in the integrable approach, one needs to trade this term for a shift 
of the field <p i n the Kondo interaction: the way this trading takes place depends on the 
charge of the exponentials, and this is why there is a rescaling in the TBA expressions, 
which are formally computed using an action with a H dependent Kondo coupling: see 
|23|1 and below for more details. Since all the integrals near the IR fixed point are defined 



by analytical continuation, they clearly lead to results obeying (|5.9|). In addition, the 
local and non local integrals of motion commute: therefore, in physical properties that 
involve the logarithm of the partition function (for instance), the terms coming from the 
Kondo type perturbation near the IR fixed point do not mix with the terms coming from 
the integer spin conserved quantities i. Therefore, to the non analytic contribution to / 
near the IR fixed point, is simply added an analytic contribution in odd powers of H/Tb- 
The structure of this analytic contribution is actually extremely simple, and depends only 
weakly on the spin. 

Of course, the argument establishing the duality also holds at non vanishing temper- 
ature. Even though no close expression is known for the free energy in that case, we thus 
expect (|5.9[) to still hold, this time up to terms analytical in powers of H/Tb-, T/Tb- 



Observe that fl5.8|) is very similar to ( 3.23 ). It would become identical if the the boundary 
sine-Gordon term came with the coupling 2Asin7r<7, which is actually the natural choice within 
the quantum group framework underlying these problems. 

9 This remarkable property is clearly visible on the logarithms of the reflection! matrix ele- 
ments; see the formula ( |3.19| ) for a completely analogous example in the context of the boundary 
sine-Gordon model. 
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5.2. Duality in the boundary sine- Gordon model 



That (|3.22|) contains only one cosine has a simple physical meaning here - the flow 



approaches the IR fixed point along a direction where there is a term in the hamiltonian 
corresponding to tunneling of electrons, but no term for tunneling of pairs, triplets etc. 

The most interesting properties to study in that context are transport properties, for 
which a non equilibrium formalsim such as Keldysh is required. Some of our conventions 
are discussed in the appendix; here, we will concentrate on the salient features only. In- 
troducing a vector potential a(t), the current is computed by I(t) = Sl I lZ , and expanded 
perturbatively near the UV or IR fixed point. Consider the UV fixed point first: there, 
the potential vector can be reabsorbed into the cosine term by a shift of the boson, so, 
restricting to constant voltage V, the current expands as a series of Coulomb gas integrals 
somewhat similar to the ones in equilibrium; the key difference however, is that the vertex 



operators V± = exp ±iy/2ivg^ are integrated on the Keldysh contour, represented in figure 
3, and that contour ordered propagators are used. 



Fig. 3: Keldysh contour. 

More specifically, a vertex operator stands at the extremity of the contour at time t, 
while other operators are integrated on the contour. The only non vanishing contributions 
are the ones which are electrically neutral. At non vanishing voltage V, each vertex 
operator V±(t') comes with an additional phase e l±gVt . The contour ordered propagator 
is such that (T c [<I> (£')$(£")]) = — ^ ln(t> — t<), where t> denotes the time that is the 
latest as measured along the contour, and t < the time that is the earliest. For instance, 
for t' on the above or lower part of the contour as in figure 3, the contraction that appears 
in the computation of the current is, (T c e W^9^W Q-W^g^it ) ^ = ^ 2g ; resp. nT^pg • 
The non trivial monodromy of the vertex operators ensures that the contribution of the 
two parts of the contour do not cancel out, and a non trivial result is obtained (see the 
appendix for some examples). At finite temperature T, the only change is that, in the 
propagator, hit is replaced by In sm ^ Tt . 

Now, our point is not so much to discuss the structure of this expansion (many details 



on this issue can be found in [23J for instance), but to comment on the duality proper- 



ties it might give rise to. For this, let us investigate the computation of the current in 
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the IR. If, in our framework of dimensional regularization, the only operator in the IR 
were the cos duality would easily follow from the matching of the two expansions. 

The complication we have to discuss is the role of all the O2/C+2 operators added to the 
hamiltonian. 

Let us first make a crucial observation. Consider for instance the contour ordered 
propagator in the two situations of figure 4, 




t" 




r 



Fig. 4: Conserved quantity on the contour. 

where say <9$(t') (we call in this paragraph $ what is in fact the dual of the original 
field $ for notational simplicity; d denotes time derivative) is inserted at time t' , and some 
expression X(t") at time t" . In the first situation where t' occurs earlier on the contour, 

(T c [dHt')X(t")]) = J2 2n(t»-t>) j[{n 

contractions 

where X denotes the remainder in X once contracted. In the second situation where t 1 
occurs later on the contour, 

(T c [d*(t')X(t»)]) = £ 

contractions 

and of course, the two expressions are actually equal. This easily generalizes to cases where 
9$ is replaced by any polynomial in derivatives of $, in particular the C>2k+2- As for X, it 
can be one of the C2/C+2 itself, as well as a product of such an operator by a vertex operator, 
the result still holds: in other words, the result of the contraction does not depend on the 
order on the contour, and for the C2/C+2 operators, integrals along the Keldysh contour 
just behave like ordinary contour integrals - a somewhat trivial fact, once one remembers 
that there is no cut in the complex plane for the contractions involved here. 

This observation being made, we observe that the 02fc+2 in the IR hamiltonian give rise 
to two complications. First, when one trades the coupling of the vector potential J adt<& 
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for a shift of the field $, this time not only does one get a shift -Vt in the argument of the 



cosine; one also gets in the new hamiltonian additional terms made of a and polynomials 
in the derivatives of $ (for instance, the : (9$) 4 : in : T 2 : gives rise to a a 2 : (9$) 2 : term, 
etc). Thus, when defining the current as the functional derivative of InZ with respect to a, 
one gets, in the IR, a more complicated expression than in the UV: what has to be inserted 
at t on the Keldysh contour is the sum of a vertex operator and a series of polynomials in 
derivatives of 

Now consider the perturbative computation of this current in the IR: we have to insert 
on the Keldysh contour either vertex operators or operators C?2fc+2- For the component of 
the current at t that is not the vertex operator however, no cut is necessary at t. According 
to the observation above, the integrals on the Keldysh contour of the various insertions 
then just behave like ordinary integrals, for which the upper and lower parts of the contour 
cancel out - in other words, the current is still obtained by only inserting vertex operators 
at t, in complete analogy with the UV case. 

The second complication due to the 02fc+2 is that these operators contribute to the 
perturbation series in the IR. Consider thus a generic term in the perturbation series, where 
a few vertex operators as well as conserved quantities have been inserted. To regulate 
divergences, it might be necessary to slightly displace the contours - this does not matter 
anyway, as we now argue. Indeed, consider moving the contours for the insertions of 
conserved quantities, say 02fc+2 and 02fc'+2- Since they are polynomials in derivatives 
of according to our observation above, these contours can be deformed as for usual 
integrals. The residue of their short distance expansion is a total derivative, so when 
we move one contour through the other, we are left with the contour integral of a total 
derivative. If in turn we try to deform this contour to zero, since the short distance 
expansion of a total derivative with any quantity cannot have a simple pole, no obstacle 
is met. In other words, we can freely pass through one another the contours for conserved 
quantities 02/c+2- 

Let us now try to pass these contours through the vertex operators. Consider thus 
a situation as the one in figure 5 where we have four vertex operators inserted on the 
Keldysh contour, and are trying to pass the 02fc+2 contour through them. 



02k+2 



t1 



t2 




t3 



Fig. 5: Vertex operators. 
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In doing so, we encouter four poles, whose residues are total derivatives. Let us call 
the residue of the expansion of 2 k+2 and V e , dOk, £ - If h, t 2 , £3, t are the arguments of the 
four vertex operators, the total quantity picked up is 



dt 1 O k , tl V e2 (t2)V e3 (t 3 )V e (t) + permutations 



Instead of contour integrating this quantity, let us simply look at its contour ordered 
average. Because the various contractions depend only on the difference of arguments, the 
effect of summing over permutations is to compute the derivatives with respect of sums of 
arguments of quantities that depend only on their differences, that is, is zero. Hence, the 
contribution of the residues when moving the 2 k+ 2 contour through the vertex insertions 
cancels out, and we can squeeze this contour to zero. In other words, the 2 k+ 2 do not 
contribute to the current at all. This is independent of the voltage or the temperature. As 
far as the current goes therefore, it is fully determined, in the scheme where integrals are 



analytically regularized, by the 2A cos perturbation in the UV and 2AdCos \J^f& 



in the IR. This allows us to conclude that 



6. Conclusions 

In conclusion, we would like to stress that the implementation of the IR perturbation 
theory, as well as the existence of duality, rely completely on the integrability of the 
problem. The latter acts as a symmetry that restricts the IR hamiltonian in a drastic 
fashion, so that the structure of the IR perturbation is almost the same as the UV one, 
maybe up to analytical terms. In general impurity problems, we do not expect the duality 
to be more than a quick qualitative argument to find the leading irrelevant operator. We 
also do not expect IR perturbation theory to make much sense, because of the difficulty 
in regularizing higher order terms when the operators do not commute. 
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Appendix A. Normalization of conserved quantities in the sine-Gordon model 



We discuss here the problem of determining the constant \2k+i in the definition fl3.4|) . 
To do so, we consider the free action in the bulk, to which we add two perturbations: one 
of them is a field coupled to the U(l) charge, and the other a term proportional to the 
conserved quantity X 2 k+i- Going to a hamiltonian description in the closed string channel, 
we obtain 

[2a r°° r°° 
n = n free + vJ^- dyd y $ + (jl / o 2k +2. (A.i) 

V " J — oo J — oo 

The question we then consider is the ground state energy of this theory. It can be computed 
using the integrable structure [p5|,|[26||,P|. We use here the notations of section VI of the 



latter reference; in addition, we set M = 2, ft = 1, V = 2V and e = 1. We denote the 
rapidity by (3 instead of 9. The constant A of this reference corresponds, in the present 
paper, to A = - — 1; it of course has nothing to do with what we call A in the present 
paper, that is the bare coupling. 

The problem factorizes into R and L components. Consider say the R sector. The 
equations determining the ground state density of particles depend only on the momenta, 
which are not affected by the perturbations of the hamiltonian. The cut-off A however 
changes, in a way that depends on the perturbations in a crucial way. The equivalent of 
eqn (6.9) of is now 

V - e 13 - ^hh±l e W+i)p = e(/3) _ r _ p) e (j3') d p'. (A.2) 

2 J -oo 



-iujA 



It follows, if the Fourier transforms are defined as in ||, and if we denote e_(a;) = we' 
that 

% u-i i u , 

liX 2k +i G-(u)G+[(2k + l)i] (2k+1)A [ ' } 

2i u-(2k + l)i 

where the kernels are given in ( |5.3|) above, and in eqn (6.6) of 0. The cutoff A is such that 
e(A) = ie lim^^oo we_ (iu) = 0. We will restrict ourselves to the case where V, fx » 1, 
where the first term in (|A.3| ) becomes negligible. It then follows that 

(2fc+l)A _ 2V G +(°) (A4) 

V\ 2k+1 G + [(2k+l)iY 1 • ] 
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We can then compute the energy per unit length 



E = 2 dPp(P) \e? + fiX 2k+1 e^ k+1 ^ -v]n 2p\ 2k+1 - 2Vp{Q)p[-{2k + (A.5) 



where again we neglected the first term, of order one. Using another result form [f 



p(u) 



1 G-(u)G+(i) 

2Z7T u — i 



(iu>+l)A 



(A.6) 



we get, after some algebra 
E = 



1 2k + 1 . _ 2fc+2 

V 2k + 1 



Ti 2k + 2 



^ G+(i)G + (0)w 
G + p + l)# 



(A.7) 



On the other hand, we can compute the energy directly from the hamiltonian. In the 
integral of 2k + 2 , only the leading term contributes since all the others involve second or 
higher derivatives of <fi, which vanish at the saddle point. Therefore one has, using the 
normalization in (|3.8| ), 



2 fc + 1 



,2q2k + l 2 k +2 /2\ 2fe + 1 
E = -2V\ — — —V^ 1 



it 2k + 2 
From this, it finally follows that 



1 



A 



2fc+l 



7T 



- (k+iy. 

9. 
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<?)J 



2k + 2 



2k+l 



(-27T) 



2 k + 1 
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(2fc+l)g~ 


2(1-3) J 
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2fc+l 


1.2(1-0) J 
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Appendix B. Some remarks on Keldysh and analytic continuation 



In ||27|1 , a formula for the current was proposed 



I = gV- ig nTj-\n Z ^ X) 



9 In A Z(-p,X) 



(B.l) 



Here, Z(p) is an analytic continuation of the partition function at "imaginary voltage", 
accomplished through a Jack polynomials expansion ||27|| . 

The partition function at imaginary voltage is defined as Tre~' H ^ p ^ T , where p is 
an integer and H.{p) is obtained from 7i by shifting the argument of the exponential 
cos ^2ng<& — > cos ( K \f2ng<fr + 2ivpTy). The physical voltage is such that 2npT = igV, so 
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an analytical continuation in p from integer to imaginary values has to be carried out. The 
Keldysh formalism actually tells us how to perform this continuation. To see this, let us 
first recall some basic results. We consider the boundary sine- Gordon model with a vector 
potential a(y) @. After the usual shift, one can write the partition function 

1 + A 2 J J d yi dy 2 cos[o(j/i) - a(y 2 )]p(\ yi - y 2 \) + . . . L (B.2) 

where dots stand for higher order terms, we restrict to g < | so the integrals are all 
convergent, and 

sm 7rTy / 

The current then follows from I = 5 l n z , together with the usual contour deformation 

I(t) = 2A 2 / — sin[a(£) - a(t')] {Tj^^^e^^^) + ... 

1 (B4) 
nl2 /■* ,,. r/ , , A1 P>(t-t') -P<(t-t') K ' } 
= 2A 2 / <ft' sin[a(t) - a(t')] — + • • • , 

J — oo * 

where P > (resp.P < ) is the analytic continuation to y = it(resp.y = —it) of P. Let us now 
restrict to a DC voltage a(£) = gVt; in that case, 

f = A! Vhi#), (B.5) 



where 



One has 



, tot P>(t)-P<(t) 



P(x) = / cfte lxt ^. (B.6) 



o 



"OO J j. 

P(x) = 2(7rT) 29 " 1 sinvr^ / e™ l ^ T — 



o 



(sinht) 



23 ' 



The latter integral is tabulated, so one gets 



P(gV) = (2 ? r) 2 ^- 1 ; in%9 . x — F(1 ~ 2 4 (B.7) 



Let us now get back to the question of analytically continuing Z(p) in (|B.1|) . Consider 
the case p is an integer, so the quantities Z(p) are well defined. Of course, for p integer, Z is 
even in p, so the argument of the derivative in ( p.lj ) is identically zero. Let us nevertheless 
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expand it formally in powers of A. At lowest order, the contribution to the current involves 
the quantities 

,1/T 

Q(p) = / cos(27rpTy)p(y)dy. (B.8) 



One has 

dy 



Q(p) = (ttT) 23 " 1 ( e 2ipy - 
Jo 



(smy) 29 

This integral is tabulated, and, for p an integer, reads 

Q(p) = (2it) 2s T 2 ^f-l)P- r( \~ 2 ^ 

VW 1 1 1 1 r(l-g+p)T(l-g-p) 



(B.9) 



Still for p an integer, the expressions ( |B.7| ) and (|B.9| ) coincide: this is because, in that case, 
the voltage plays the role of a Matsubara frequency, so what we have here is the standard 
identity between Fourier transforms of temperature Green functions and retarded time- 



Green functions [28]. But we see now how to perform the continuation of Z(p): we first 
need to deform the contour, at each order in peturbation theory, from the imaginary time 
interval to the Keldysh contour - this is possible for p an integer - and then replace p by 
igVT/2n in the integrand. 

At the level of final expressions, that is ( ]B.7|) and ( ]B.9|) , however, it is less clear what 
must be done. The "recipe" proposed in |2JJ consists in expanding the integral Q into a 
sum of rational functions of Gamma functions 

and then perform the continuation simply by replacing p by the appropriate non integer 
value in each of the Gamma functions. The sum ( |B.1U| ) was studied in details in [ |27[ , where 
it was shown that, for arbitrary p, this continuation of Q coincides with ( p.7| ). Therefore, 
this gives the same result as the one obtained by deforming contour, which is the correct 
physical one, based on the Keldysh analysis. 

Notice that the continuation using the expression ( [B.IO] ) is not the same as the con- 
tinuation discussed in ||29|| . In the latter work, the author expands the partition functions 
Z(p) in (|B.1|) over Matsubara propagators, and then performs the continuation in p in each 
term independently. There is no reason why this definition should coincide with ours, and 
herefore it is not surprising that disagreements are found in J29| : only one continuation is 
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expected to work, and the conjecture made in [57] is that, for the tunneling problem de- 
scribed by the boundary sine-Gordon model, it is the one using the expansions ( p.lOj ), and 
for higher order, the corresponding sums based on Jack polynomials theory (see below). 

To understand better why this might be true, we observe that P(gV) has a simple 
power law behaviour P oc (V/T) 29 ~ l as V/T — > oo, while Q(p) in ( p3 . 9|) does not. The 
power law behaviour is expected from the Keldysh contour representation, and on physical 
grounds as well: it is necessary for the current to have a finite expansion in terms of V/Tb 
as T — > 0. On the contrary ( [B.lOp , supplemented by the naive replacement of p by 
does have the right behaviour. 

When one considers higher powers of A in the expansion of the current using the 
Keldysh formalism, one gets integrals which still coincide, for p integer, with the integrals 
occuring in Z(p), through contour deformation. The challenge, if one whishes to prove the 
conjecture (|B.1| ), is to show that replacing p by igVT/2n in the deformed contour integrals 
coincides with replacing p by igVT/2n in the expansion 

Q 2n ( P ) = ( 27r) 2 gr2 n g - 2 n + i i y tt r [m + g(n - i + i)j 

K) 1 ' r(g)^^l\r[m i + l+g(n-i)} 

r [p + rrti + g(n — i + 1)] 
r \p + rrii + 1 + g(n - i)] ' 



x 



(where the sum is over all sets m = (mi,...,m n ) with mi > m<2 ■ ■ ■ > > 0), a 
result we explicitely checked above at lowest order. Clearly, the two procedures define 
analytical continuations of functions defined for integers. Under reasonable assumptions, 
it is known that two functions that coincide on integers and have the same behaviour 
at infinity are actually identical. Therefore, a proof of ( P . 1|) would simply be that the 
continuation of ( |B.11|) behaves, when p — > ±zoo, as p 2n 9- 2n + 1 . Though we have not 
proven this analytically, we have checked it numerically for the first few values of n. 
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